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Monte Carlo simulations have been used to study magnetic ordering in coupled anisotropic 
ferro/antiferromagnetic (FM/AFM) films of classical Heisenberg spins. We consider films with 
flat interfaces that are fully uncompensated as well as rough interfaces that are compensated on 
average. For both types of interfaces above the "Neel temperature" we observed order in the AFM 
with the AFM spins aligning collinearly with the FM moments. In the case of rough interfaces 
there is a transition from collinear to perpendicular alignment of the FM and AFM spins at a lower 
temperature. 



Ferromagnetic/antiferromagnetic (FM/AFM) bilayers 
have been shown to exhibit a unidirectional shift in the 
hysteresis loop (exchange bias) and a significant increase 
of the coercivity 0. Magnetic properties of FM/AFM 
bilayers can also be drastically different from those of 
free FM and AFM films. For example, even though in 
most exchange bias systems the blocking temperature is 
lower than the Neel temperature of the antiferromagnet, 
order in the AFM has been observed above the Neel tem- 
perature due to the coupling to the ferromagnet Q. 

Earlier theoretical work 0, 0] has shown that for com- 
pensated interfacial exchange (i.e. zero net exchange in- 
teraction across the FM/AFM interface) the FM aligns 
perpendicular to the AFM easy axis. Although similar 
perpendicular orientation has been observed in numerous 
FM/AFM systems 0, understanding the nature of the 
interfacial coupling and roughness in FM/AFM bilayers 
remains a challenge. 

In this paper we use Monte Carlo simulations to study 
the effect of interfacial exchange interaction and rough- 
ness on transitions from ordered to disordered states in 
FM/AFM bilayers. The work is motivated by recent ex- 
periments of Ijiri et. al. [jg, which show for CoO/Fe3C>4 
multilayers, that the transition to perpendicular ordering 
takes place at temperatures considerably lower than the 
FM and AFM ordering temperatures. 

The model studied here consists of a ferromagnetic 
(FM) film coupled to an underlying antiferromagnetic 
(AFM) him where the lattice is coherent across the 
FM/AFM interface. The structure of the Alms is a body- 
centered cubic (BCC) lattice, with linear sizes L x = L y = 
L < 96. Each him is composed of 12 staggered (because 
of the BCC structure) layers of classical unit length spins 
S r = (S* , S% , S"), which interact via the Hamiltonian 
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where (r, r') denotes nearest-neighbor pairs of spins cou- 
pled with exchange interactions Jf > in the FM him, 
J A < in the AFM film, and Jj at the FM/AFM inter- 
face. Spins in the AFM film have a uniaxial single-site 
anisotropy Ka > 0, whose easy axis is along the y axis. 
The demagnetizing field on the FM film is modeled with 
a hard-axis (Kf < 0) along the z direction, which is per- 
pendicular to the FM/AFM interfacial plane. No exter- 
nal magnetic field is applied. We use periodic boundary 
conditions along the x and y directions and free boundary 
conditions along the z direction. 

We consider interaction and anisotropy parameters 
Jf = 5J > 0, Ja = —J, Ji = —J, Ka — J, and 
Kp = —0.5 J. We model flat interfaces as well as "rough" 
ones which are comprised of uniformly spaced steps, with 
6, L, and one spin per terrace in the x, y, and z directions, 
respectively. The terrace sizes along the x and z direc- 
tions are kept fixed for the different him cross sections 
used. For both types of interfaces the exchange coupling 
across the interface is uniform (Ji = — J); therefore, the 
flat interfaces are fully uncompensated while the rough 
ones are compensated on average. 

The FM and AFM order parameters are the uni- 
form and staggered magnetization per spin, respectively. 
These are defined as m = | S r | /Nf and m s = 
I Sre/ ^ — J^reii ^r\/^A, where I and II denote the 
two simple cubic sublattices of the BCC lattice and Np 
and Na denote the number of spins on the FM and on 
the AFM films, respectively. The summations here are 
performed over all sites on both types of films. The mag- 
nitude of the two components of the uniform magnetiza- 
tion parallel to the interfacial plane is denoted as m x and 
m y . Our simulations were carried out using importance 
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sampling Monte Carlo methods , with Metropolis algo- 
rithm, at fixed temperature T. Typically 2 x 10 5 Monte 
Carlo Steps/site (MCS) were used for computing aver- 
ages after about 3 x 10 5 MCS were first discarded for 
thermalization. Whenever not shown, error bars in the 
figures are smaller than the symbol sizes. 

FigG] shows the uniform and staggered magnetization 
as a function of T for systems with a rough (stepped) 
interface. At low T, both m and m s tend to unity, in- 
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FIG. 1: Uniform and staggered magnetization as a function 
of temperature for L — 24, 48, and 96, with rough interface 
and parameters Jf = 5 J > 0, J a = Ji = — J , Ka = J, and 
K F = -0.5J. 

dicating ordered FM and AFM films. At very high T, 
both m and m s decay to zero, indicating disordered spin 
configurations on both films. Finite-size effects are illus- 
trated using films with cross sections L = 24, 48, and 96, 
where the number of spins per terrace along the x di- 
rection is constant. These data suggest the existence of 
a FM transition between ordered and disordered states 
at T c w 9.3J/ks and a possible Neel temperature at 
Tjy ~ 2.2J/ks, which coincides with the Neel temper- 
ature of a free 12-atomic-layer AFM film 

Temperature dependences of m x and m v are shown in 
Figs|2^ and \Bp for systems with a rough and a flat in- 
terface, respectively. For both types of interfaces, the 
spins on the FM film orient predominantly in a direction 
collinear with the easy axis of the AFM at high temper- 
atures, even above TV. This is an indication that there 
is still order in the AFM above Tjv due to the coupling 
to the ferromagnet. As the temperature is lowered in the 
case of the rough interface the FM spins switch to ori- 
ent in a direction that is perpendicular to the AFM easy 
axis, which is the direction of the AFM spins. Our results 
suggest that the onset of this perpendicular orientation 
is very sharp and it occurs at a temperature below the 
Neel temperature. The ^-components of the uniform and 
staggered magnetization are very small for all T, Note 
that in the absence of the AFM film, spins on the FM film 
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(a) L=96, rough interface 
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FIG. 2: Temperature dependence of the uniform and stag- 
gered magnetization and of the magnitude of the components 
of the uniform magnetization m x and m y for (a) L = 96, with 
rough interface and for (b) L = 48, with flat interface. Both 
cases are for Jf — 5J > 0, Ja = Ji = —J, Ka = J, and 
K F = -0.5J. 



have global rotation symmetry in the x-y plane, which is 
the easy plane for spins on this film. The preferential ori- 
entations of the FM spins observed here either below or 
above Tjy result from the exchange coupling to the AFM 
film. 

In studying thin films, it is also interesting to examine 
the layer magnetization since this can vary significantly 
depending on the proximity to the interface or to the free 
boundaries. Figs|3^, and \S]p show the layer magnetiza- 
tion and the magnitude of its y-component, respectively, 
as a function of the film layer for a system with rough 
interface and L = 96. In our notation, layers 1 to 12 
comprise the AFM film, layer 13 has uniform terraces 
belonging alternately to the AFM and FM films, and 
layers 14 to 24 form the FM film. For each temperature, 
the slightly lower magnetization on the surface layers is a 
consequence of the free boundary condition along the z- 
axis. In addition, the compensated exchange on average 
across the interfacial plane reduces the layer magnetiza- 
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FIG. 3: (a) Layer magnetization and (b) the magnitude of 
its y-component as a function of the film layer, for rough 
interface and L = 96. Parameters used are Jf — 5 J > 0, 
J a — Ji — —J, Ka — J, and Kf = — 0.5J. 



tion for the layers near the interface. The perpendicular 
orientation of the FM spins with the AFM spins at low 
T can also be seen in Fig|3jD, where e.g. for T = 0.6J/fcs 
the magnetization on the FM layers has a very small y- 
component. At a higher T, e.g. T = 1.8J/kg < 2jv the 
magnetization on the FM layers has a large y-component. 
When the interface is flat with uncompensated exchange 
across it, the reduction of the magnetization on the layers 
near the interface is not very large, at low T (see Fig0J 
where layers 1 to 12 comprise the AFM film and layers 
13 to 24 the FM film). 

Simulations using Jp = 5 J > 0, J a = — J, Ji = — 2 J, 
Ka — J, and Kp — —0.5 J {i.e. with a stronger in- 
terfacial coupling) with rough interface suggest that the 
onset of the perpendicular orientation of the FM with 
the AFM easy axis occurs at a lower T when the interfa- 
cial coupling increases in magnitude. We also find that 
for a film with a rough interface that is compensated 
on average, when the interaction and anisotropy param- 
eters are such that the AFM becomes disordered at a 
higher temperature than the FM (e.g. with Jp = J, 
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FIG. 4: Layer magnetization as a function of the film layer, for 
flat interface and L — 48. Parameters used are Jf = 5J > 0, 
Ja = Ji = -J, K A = J, and K F = -0.5 J. 



J A = —J, Ji = —J, Kp = — 0.5J, and Ka — J), the 
perpendicular orientation of the FM and AFM spins is 
present throughout the FM phase. 

Besides the usual critical slowing down near critical 
points, studies of phase transitions in thin films are fur- 
ther complicated by the fact that order in the films varies 
layer by layer, as illustrated above. Therefore, determin- 
ing the nature of phase transitions may require a layer 
order parameter. However, the goal of our current work 
is not to study the thermodynamic limit of these bilay- 
ers, but rather to understand experimental observations 
of finite systems. 

Monte Carlo simulations have been used to study a 
system of coupled FM/AFM films, with flat interfaces 
that are fully uncompensated as well as rough interfaces 
that are compensated on average. For both types of in- 
terfaces above the Neel temperature we observed order 
in the AFM with the AFM spins aligning collinearly with 
the FM moments. As the temperature is lowered in the 
case of the rough interface, we have seen a transition from 
collinear to perpendicular alignment of the FM and AFM 
spins. 
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